Self-localization of holes in a lightly doped Mott insulator 
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We show that hghtly doped holes will be self-trapped in an antiferromagnetic spin background at 
low-temperatures, resulting in a spontaneous translational symmetry breaking. The underlying Mott 
\l ' physics is responsible for such novel self-localization of charge carriers. Interesting transport and di- 

electric properties are found as the consequences, including large doping-dependent thermopower and 
dielectric constant, low-temperature variable- range- hopping resistivity, as well as high-temperature 
strange-metal-like resistivity, which are consistent with experimental measurements in the high-Tc 
' cuprates. Disorder and impurities only play a minor and assistant role here. 
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I. INTRODUCTION 



There exist several types of electron localizations in condensed matter physics. In the presence of disorder or 
impurities, waves can be localized due to quantum interference, which is generally known as the Anderson localization 
^ ' 0. In a two-dimensional (2D) system carriers are expected 0] to be always localized with the resistance diverging 
"j^ I either logarithmically ("weak localization") or exponentially ("strong localization") as T ^ 0. A different kind of 
i: localization involves the self-trapping of small polarons in strong electron-phonon interacting systems. Both types of 

localizations here mainly concern non-interacting or weakly correlated electrons. 
, In the high-Tc cuprates, the undoped system is a Mott insulator, in which the charge degree of freedom is totally 
^ frozen out by strong on-site Coulomb repulsion. Its spins form an antiferromagnetic long range order (AFLRO) at 
O , low temperatures, which quickly collapses upon hole-doping .3|. At small concentration, S < 0.05, the doped holes as 
I ^ I ■ charge carriers remain localized with the low-T resistivity well fit 0, 0, IS by those of traditional variable-range- 
' hopping (VRH) type, usually applicable to doped semiconductors. Only when the hole concentration exceeds 0.05 
• will the charge carriers be truly delocalized, where the ground state becomes a d-wave superconducting state. 
' But lightly doped high-Tc cuprates are of no conventional doped semiconductors. Here they are doped Mott 
: insulators with the majority of the charge degree of freedom still remaining frozen. The strongly correlated effect 
should thus play a crucial role in the charge transport. The issue why doped holes should be always localized at small 
J doping and how they eventually become delocalized with increasing doping is a very important question, which may 
' also be relevant to understanding the microscopic mechanism of high-temperature superconductivity in the cuprates. 
Tj" I Recently, a new kind of self-localization, which can be particularly attributed to the Mott physics, has been proposed 
■ Hi ^3 J in which the pure Coulomb interaction is responsible for the self-localization of electrons near half- filling. 
Here disorder or impurities is no longer as essential as in the Anderson localization, although they may be helpful to 
_i . induce a true translational symmetry breaking at T = 0. Such a self-localization of holes is caused by the phase-string 
^ ' effect referring to a stringlike defect left by the hopping of a hole, which is irreparable in the ground state of 
the t — J model no matter whether there is an AFLRO or not at arbitrary doping. Furthermore, it has been also 
^ ' shown that the charge localization does not contradict to the photoemission experiments |0 in which the observed 
Q ] "quasiparticle" dispersion in the single-particle spectral function can be well accounted for in terms of the "spinon 
O ■ dispersion" . 

^ ■ In principle, it is not surprising that due to the separation jl3| of spin and charge degrees of freedom in a doped 
Mott insulator, the charge carriers will get localized at low doping by scattering with fluctuations from an independent 
of degrees of freedom. In fact, in a different gauge-theory approach to the t — J model, the localization of charge 
,, carriers has been also obtained li"!, by scattering to the gauge fluctuations. 

In this paper, we study such a self-localization phenomenon based on an effective phase-string description |llL Il5| | 
of the t ~ J model. We shall focus on the spin ordered phase (not necessarily long-range ordered) at very low doping 
and explore the unique charge self-localization behavior. As shown previously, a doped hole in such a regime will 
induce ^ J^] a dipole-like spin structure. We find that the kinetic energy of the hole is severely frustrated by the 
phase string effect such that the hole-dipole is self-trapped in real space. We discuss the corresponding transport 
properties and show that the thermopower saturates to a Heikes-like formula as the result. While the resistivity 
exhibits a Mott-VRH-like behavior at low temperatures, the collapse of the hole-dipole composites and the release of 
free "holons" at high-temperatures will lead to a strange-metal-like behavior there. Furthermore, the existence of the 
hole dipolar structure predicts a large doping-dependent dielectric constant, which diverges at the deconfinement. All 
of these properties seem to paint a consistent picture for the complex transport and dielectric phenomena observed 
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in the high-Tc cuprates. 

The remainder of the paper is organized as follows. In Sec. II, we discuss the self-localization of doped holes in the 
spin ordered phase at low doping based on the phase-string model. A renormalization group (RG) analysis will be used 
to determine the phase diagram. In Sec. Ill, we discuss the experimental implications of the self-localization, including 
the thermopower, resistivity and dielectric constant, and make comparisons with experimental measurements. Finally, 
the conclusions are presented in Sec. IV. 

II. SELF-LOCALIZATION OF HOLES IN A DOPED MOTT INSULATOR 

We shall adopt the phase-string model as the microscopic description of how doped holes move in an antiferromag- 
netic (AF) Mott insulator. This model is obtained ^5] as a low-energy effective theory based on the t — J model, 
which can accurately describe AF correlations at half-filling and possesses a d-wave superconducting ground state 
at doping concentrations larger than Xc (— 0.043 at T = 0) 0. What we will be interested in the following is 
the non-superconducting phase below the critical doping x^, where an AFLRO or a spin glass state persists at low 
temperatures. 

The existence of such a non-superconducting phase is the consequence that the long-range AF correlations (not 
necessarily AFLRO) win in the competition with the kinetic energy of holes at sufficiently low doping. In the phase- 
string model, this phase will be characterized by the "spinon condensation" . In the following, we shall analyze in detail 
the behavior of doped holes in this low-doping regime and demonstrate that the charge carriers must be self-localized 
in the ground state, resulting a spontaincous translational symmetry breaking. 

A. Phase-string model 

We start with the phase-string model -ffstring — Hh + i?s, which is composed of two terms: The charge degree 
of freedom as characterized by the "holon" term 

where th ^ t and the "holon" operator, hj, is bosonic; The spin degrees of freedom as described by the "spinon" term 

iJ, = -J,^ (e-<)6,L&t_^+iJ.c. (2) 

where Js ~ J and the "spinon" operator, foj^, is also bosonic. 

Basic features of this model are as follows. At half filling, the gauge field A'lj can be set to zero in Eq.® and 
reduces to the Schwinger-boson mean-field Hamiltonian which describes both the long-range and short- 
range AF correlations fairly well. Upon doping, A'^j is no longer trivial as it satisfies a topological constraint: 
J2c "^ij ~ ""SiGEc (^C" denotes the area enclosed by C) with denoting the "holon" number at site I, which is 
interpreted as that each "holon" behaves like a 7r-fluxoid as felt by the "spinous" . Thus, Af^ will play the key role of 
frustrations introduced by holes that act on the spin degrees of freedom. Similarly, the "holons" are also subjected 
to frustrations from the spin background, through the gauge field Afj in Eq.lQJ. Here Afj satisfies a topological 

constraint: X^c ~ ""X^jeSc (""^'t ~ ""u) ^^^^ ^L- denoting the "spinon" (with index cr) number at site /, which 
can be interpreted as that each "spinon" behaves like a ±7r-fluxoid as perceived by the "holons" . 

The spin and charge degrees of freedom are thus mutually "entangled" in the phase-string model i/string based on 
two topological gauge fields, A'^^ and A^j. It has been shown that if the holons (which are bosons) experience a 
Bose condensation at larger doping, the spinous, which behave as vortices according to Afj, must be "confined" at 
low temperatures. Such a "spinon confining phase" actually corresponds to the d-wave superconducting phase in the 
model. Here, the spinon confinement occurring in the spin degree of freedom is closely related to the superconducting 
phase coherence in the charge degree of freedom. 

In contrast, if the spinous (which are also bosons) are Bose condensed at low doping, free holons as "vortices" 
cannot live alone either, and must be also "confined" 'o', 'lO]. Such a spinon condensed phase is thus also known as 
"holon confining phase", which constitutes a non-superconducting solution of the phase-string model in the lightly 
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FIG. 1: Holes are self-trapped in real space a,t S < Xc by the phase-string effect, carrying dipolar spin configurations. Holes 
become delocalized at 5 > Xc, with the critical doping concentration Xc — 0.043 (see text). 

doped regime. In the following, we will demonstrate how holons will get self-trapped in real space in such a phase. 
The spin ordered phase as characterized by the spinon condensation is therefore intimately related to the charge 
self- localization (holon confining) . Fig. 1 schematically illustrates the above-mentioned two phases at different doping 
concentrations . 



B. Antimerons induced by holons: "Flux-quantization" condition 

Let us recall that at half-filling, the Bose condensation of spinous, 

< >^ (3) 

will naturally gives rise to an AFLRO lying in the x-y plane with (S^) = (—1)* < >< bti >. However, once a 
hole is added into such an ordered state, the energy cost associated with a bare holon, if the condensation (PJ remains 
unchanged, would become logarithmic ally divergent in terms of 

^ I J,ln(L/a) ^ oo, (4) 

where Jg = Js < bi^ >< b^ >, A'' is the continuum version of A'^ , L is the size of the sample and a is the lattice 
constant. 

But this is not the correct physical picture. In reality, the two-component spinon "superfluid" in Eq.Q can response 
to the presence of holons by forming a spin supercurrent. Such a spin supercurrent can then screen out the effect of 
the "magnetic fluxes" introduced by A'^j in Hg, just like the flux-quantization phenomenon in a superconductor, such 
that the resulting renormalized holon will acquire a flnite self-energy. In the following, we show how this screening 
effect takes place in the phase-string model (0. 

Define < bi^ p^a'^Zia-, ifi G A sublattice, and < bia >= y^p^a'^z*_^, if i G i? sublattice, with |zi|p-f = 1. 
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Here denotes the spinon "superfluid" density. Then the condensed part of Hg can be written as 

^cond ^ ^ ^^^.A'l,^^ < ^t^ ^t_^ > ^^^^^ 

^ -^^ E E (i + '■'^^ - \ (4)') + ^-^^ 

~ S° + iJs j Sv^aA^ ■ iz*Vz^ - c.c.) + Js J d^r (A^)^ (5) 

in the continuum Umit, where Js = JsP^o-^, E^! = —AJsN, and NN denotes the nearest neighboring sites. The spin 
supercurrent, defined by J" = —dHf^^^'^/dA^, then reads 



J" = 2 J, (v" - A'') (6) 

where 

v^^-^^a(z;Vz.-c.c.). (7) 

(T 

Generally one expects a complicated distribution of the spin current nearby a holon. But in a distance far away 
from the holon, the supercurrent J" should vanish in order to ensure the finiteness of the energy cost introduced by 
a holon. The same requirement has been used in a superconductor to realize the flux quantization. By contrast, here 
we are dealing with the "flux-quantization" in a two-component superfluid problem with an internal gauge freedom. 

By requiring the spin supercurrent J'' vanish at the boundary Coo such that 

dr • J" ^ 0, (8) 
one arrives at 



dr-w" - (t dr • A*" = 0. (9) 

C*co Coo 

Note that §^ dr • v" = / (Pr (V x v**) • z. By introducing a unit vector 

n =zaz 

where a is the Pauli matrix and z = (z*, z^), one can straightforwardly show that 



X v') ■ z = ^n- d^nx dyii. (10) 

Thus, Eq.|(Hl finally reduces to the following "flux-quantization" condition 

Q = J ^{n-d^n xdyu)^^ £ dr ■ A'' 

= -N'' (11) 
2 ^ ^ 

where is the total number of doped holes. 

Therefore, in the spinon condensed state, a holon will always induce a "screening" response from the spinon 
condensate, which is of topological nature satisfying the "flux-quantization" condition (|11|) . meaning that each holon 
will "nucleate" a spin "meron" configuration with a topological charge Q = 1/2. Such a meron configuration may be 
pictured as a spin vortex with the unit vector n lying in a spin x-y plane at a distance far away from the core, while. 
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near the core, the unit vector n starts to tilt away from the x-y plane and points towards the z-axis at the core center, 
which covers one half of the unit sphere spanned by n once, in contrast to a skyrmion which covers the whole unit 
sphere exactly once with Q = 1. 

Finally we note that such holon-induced merons are called antimerons in the earlier approach [Tol | because a 
holon itself also carries a meron (vortex) in the original spin space. To see this, let us recall that in the phase-string 
model, the spin operators are expressed in terms of spinon operators by \1\ 

5+ = (-l)^6,Vaexp[i$?] (12) 

and S~ = (S'/')^ , Sf — crbl^bi^. In the spinon condensed phase, one has 

(5+) = (-1)' < X 6a > exp , (13) 

which is twisted away from an AFLRO lying in the x-y plane by the vortices centered at holons as determined by 
Here is defined by 

$^ = ^Imln(z,-zz)nf . (14) 

But Ea. H13|l describes an bare holon effect which would result in a divergent self-energy of the holon as shown in 
Eq.Q. In order to compensate such a vortex configuration associated with a bare holon, the condensed spinon fields 
have to be twisted into 



(15) 

with an antimeron configuration 'di which is characterized just by Zicr in the present approach, satisfying the 

"flux-quantization" condition Therefore, in the phase-string model, a renormalized holon is a composite with 

a bare holon bound to an induced antimeron as demonstrated above , w' hile in the original spin space, it is an object 
composed of a meron (holon) and an antimeron, forming a dipole |9l.ll0j|. Two descriptions are equivalent. 



< h<y >-^< 6i<T > exp 
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C. Self-localization of holes: RG analysis 



According to the above analysis, the infinite self-energy of a bare holon in the spinon condensed phase can be removed 
by "nucleating" a topological spin antimeron configuration. This induced antimeron is an infinite-size semiclassical 
object which cannot propagate based on Hg. Thus, holons will be self-trapped around the cores of these antimerons. 
Due to the translational symmetry, these induced antimerons can be located anywhere in space and therefore will 
result in a spontaneous translational symmetry breaking in the spinon condensed phase. 

Physically, such a self-localization of charge carriers in the low-doping regime can be attributed to the irreparable 
phase-string effect created by the motion of holes, as discussed in Ref.Mil|. The phase-string model provides a 
mathematical framework to conveniently handle this effect. In this description, the locations of the antimeron and 
the holon inside a dipole constitute the starting and ending points of the motion of a holon, and the phase-string 
connecting such two points has relaxed into a dipolar picture, with a remaining branch-cut connecting two poles |lClj | . 

If one tries to move away the bare holon from the core of the induced antimeron, the uncompensated spin supercur- 
rent will increase the self-energy, representing an attractive potential which binds the bare holon to the antimeron. 
According to Refs.H[l3, such a potential can be estimated by 

y ~ g2 M /^g) 

a 

at |r| >a, where — irJ^, and r is the distance between the holon and antimeron. At |r| — > oo, V diverges in 
consistency with Eq.Q. 

When there are many holon-antimeron dipoles, one expects to see a screening effect on the confining potential by 
reducing V to VeS = -^V, where k denotes the dielectric constant. Previously it has been shown 9] that with the 
increase of doping concentration, eventually a transition at T = takes place, as k — > oo or Vcs ~> 0, from the holon 
self-localization (confining) phase to a delocalization (deconfining) phase a,t 6 = Xc — 0.043 (see Fig. 1). On the 
other hand, with the increase of temperature, neutral vortex-antivortex pairs like those in the XY model can also 
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be thermally excited, leading to a conventional contribution to the screening effect. To distinguish such two kinds of 
vortex-antivortex pairs, we shall call a dipole associated with a holon as a charged pair and the other kind as a neutral 
pair. In the following, we shall treat the screening effect due to both charged and neutral dipoles on an equal footing 
based on an RG treatment ■ contrast to the conventional Kosterlitz-Thouless (KT) theory , here the 

number of vortex-antivortex pairs remains finite even at T = 0, where it is equal to 6N (i.e., the number of holons). 

The probability for the neutral dipoles with two poles separated by a distance r is controlled by the neutral pair 
fugacity y^ir). In the conventional KT theory l23|, the initial ?/^(a) = e~^^" , (where (3 = and the core energy 
iJ" ~ Q^). But since the charged dipole number will be always fixed at 5 per site, the probability for the charged 
dipoles is no longer governed by the fugacity y^(a) = e^^^<: , (where E^} denotes its core ene rgy) . Instead, the initial 
fugacity must be adjusted accordingly and will turn out to be no longer so important [ol Il9l| . 

In an RG procedure, the contributions from the dipoles with the sizes between r and r + dr will be integrated out, 
starting from r — a. The renormalization effect is then represented by three renormalized quantities, X{r) = , 

y"^ (r) , and (r) , which satisfy the following famous recursion relations j^^, 0| 



dyh/dl = {2-^)yh, (17) 
dy,Jdl = (2-^)2/„, (18) 
dX/dl = 47r\yl+yl), (19) 



where r = aeK 

Define ¥"^(1) = yl{l) + yl{l), with Y^^ = yl{l = 0) + yl{l = 0). From Eqs.djl-Illl), we find 



^ = + —(X -X^)--^ In — (20) 



where Xq = X{1 ^0) ^ ^ (with k{1 = 0) = 1). The RG flow is then obtained from Eq.lO by 



X 



dX' 



2 

X, 



Att^Y^ + A{X' - Xa) - 27rln(XVXo) ' 



(21) 



The neutral pair fugacity can be determined by y^il) — e^^ ■fo^^"^^'^'- , which will show a similar behavior as in the 
conventional KT theory. 

What makes the present approach different from the conventional KT theory is the presence of a finite density of 

the charged dipoles (holon-antimeron pairs). Here, by noting ^'^4 d^r as the areal density of charged pairs of sizes 
between r and r + dr , we have the following constraint 

S/a^= /°°dr2^r^, (22) 

or 

S/a^ = r dr 2nr ^ - H dle-^' yl 



27r2a2 7o dl 27r2a2 

After an integral transformation, such a constraint can be rewritten as 

"^ T^kBT l_ Z""" j,„_2i „,2 „ o / „-2i 

g2 2'K'^d 



dle-''^-^l dle-^'yl (23) 



2^2^+^^ + -^/ dle~^'yl = 2 I e'^' X dl. (24) 



For f3q^ ^ 1, j/^ < yf-^{Q) ^ exp(— /Jg^) ^ q, the probability for neutral pairs remains a small number and can be 
neglected. 

The RG flow diagram of Eqs. (|17|) - (|19|l is as follows: The two basins of attraction are separated by the initial values 
which flow to X* ^ and yj^ — > 0, j/,* ^ in the limit ^ ^ 00. In terms of Eq.((2U, the separatrix of the RG flows 
is given by 

^^Jxo 4(X'-f)-27rln(2X'A)- ^^^^ 
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FIG. 2: The low-doping phase diagram determined by the RG analysis. The boundary set by T^aiS) separates the holon 
confining and deconfining phases. The Neel temperature Tn is determined by introducing an interlayer coupling. 



"Deconfining" temperature Tde. — Based on the RG equations of (|24|l and (|25|l . one can determine the critical hole 
density 5do = ^Ac^) or temperature Tdc = 7dc((5) at which the charged dipoles collapse and holons are "deconfined" 
from the bound state with the antimerons and become delocalized. 

We first consider the case at T = 0. Since neutral vortices do not exist at T — > 0, Eq. 124|l reduces to 

/•oo 

27r2,5dc = 2 / e"2i x dl ^ e"^' dX . (26) 

Jo Jo 

The critical doping 6dc at T = is then numerically determined by (5dc(0) — 0.043, which was previously obtained in 
Ref.lij as denoted by Xc- Now consider the limit at 6 ^ 0, where we approximately have X{1) ~ Xq — X{1 oo) = ^, 
and Eq. (|24|l becomes 

2^,^^2vrMo^. (27) 
2 

which gives rise to T^dS = 0) ~ ^ and Tdc{S) ~ Tdc(O) - Snq^/kB at 5 < 1. 

The holon "confining" and "deconfining" phases are separated by Tdc ((5) in the T — 5 phase diagram as shown 
in Fig. 2. In the "confining" phase, holons are self-trapped by binding to the induced antimerons. In the original 
spin space, each renormalized holon can be regarded as a dipolar object. In this regime, AFLRO or short-ranged 
AF ordering (spin glass) can still form, if the interlayer coupling is introduced: The Neel transition temperature is 
obtained as T^iS) « — j^^^^^^ —SSJ/ks, with a ~ 10~^ representing the ratio of the interlayer coupling J_l/J^, which 
is schematically illustrated in Fig. 2 by the dashed curve. Previously we have also shown that a superconducting phase 
will set in beyond S > Xc, where delocalized bosonic holons will experience a Bose condensation at low temperatures. 



III. EXPERIMENTAL IMPLICATIONS 



A. Thermopower 



The spontaneous translational symmetry breaking in a lightly doped Mott insulator has very unique experimental 
consequences. In such a system, since each doped hole can be self-trapped anywhere in space, due to the translational 



8 




FIG. 3: The thermopower determined by Ea. l|29|l as a function of doping concentration (solid curve). Experimental data are 
fromRefs. ^^Q^Q^ (see text). 



symmetry of the original Hamiltonian, there will be a large and doping- dependent entropy associated with the ener- 
getically degenerate configurations of holes in real-space distribution. Such an anomalous entropy associated with the 
charge carriers can be directly probed in a thermopower measurement. 

It has been previously known |22l l2.l|| that for charge carriers in a narrow band, when the temperature is raised to 
exceed the bandwidth, the thermopower will be saturated to a T-independent constant, which is entirely determined 
by the entropy change per added carrier. Namely, the thermopower Se at high-T will reduce to 



eT 
1 / dS 



e \dNe 



(28) 

E 



where /i denotes the chemical potential, S the entropy, and N^, the total electron number. The subscript E here means 
that the partial derivative is under a constant energy. 

The self-localization of doped holes means a vanishing bandwidth of the charge carriers and a spontaneous transla- 
tional symmetry breaking with a lot of degenerate real-space configurations in the distribution of holes. The formula 
(|28|l can be directly applied in such a case based on a general thermodynamic consideration. 

Putting Nfi hole on N lattice sites with no-double-occupancy would give rise to an entropy S — 
ks \nN\/Nhl {N - NhY; which then leads to the Heikes formula El = In ^i^, by noting = N - Nh 
in Eg. 1)28(1 and using the Sterlings approximation at ^ oo with Nh = N6. But since each doped hole is actually a 
dipolar composite composed of a meron and an antimeron located at two poles, which do not coincide with each other 
as wall as with other hole dipoles, the total entropy should be reduced by this fact. One can determine the entropy 
by assuming that two poles of each hole dipole are loosely bound such that it becomes a problem with a total number 
of 2Nh poles, instead of Nh, being put on N lattice sites without double occupancy. Correspondingly one obtains the 
following modified Heikes formula 

^r- = ^lnii-^. (29) 
e 2o 



This formula has no other fitting parameters and is a universal function of the doping concentration S. 
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The Heikes-type formula H29|l is plotted in Fig. 3, as the solid curve, together with the experimental data ob- 
tained at room temperatures in the Sr and Ba doped La2Cu04_y compounds (2^] (full square and bigger circle), 
Bi2Sr2Cai_j;Yj;Cu208+y system (cross), and from the recent measurement by Wang and Ong 26] in La2Sr 
2;Cu04_y (small full circle). Fig. 3 shows that Ea. (|29|l agrees qualitatively and quantitatively well with the experi- 
mental measurements in the insulating regime without any fitting parameters, where the experimental thermopower 
is sharply reduced from 300 /iV/K near half-filling to around in the metallic regime. The agreement of the theory 
and experiment quickly deviates in the metallic regime, where the experimental thermopower remains within a narrow 
range of ~ ±10 /xV/K in the optimal and overdoped regimes of the metallic phase and with a much prominent temper- 
ature dependence. We also caution that in Fig. 3, the doping concentrations in the data for Bi2Sr2Cai_a;Ya;Cu208+y 
[25I I were indirectly determined by the method involving the Hall effect which may be not as reliable as the hole 
density obtained in LSCO and LBCO compounds. 

Note that some modified Heikes formulae have been used [23. |25| phenomenologically to fit the magnitude and 
doping-dependence of the experimental data in the same low doping regime. However, it has long been a puzzling 
question why the hole bandwidth should be shrunk to an order of magnitude smaller than the temperature scale 
^ 100 K in order to explain the experiment. The self- localization of doped holes in the present theory, on the other 
hand, naturally explains this. 

An another important experimental fact is that Se has been generally found [23. 125I |2^ to decrease continuously 
to zero as T is reduced below 100 K . Such a phenomenon can be easily understood in our theory as follows. Since 
the holes are self-trapped in space, any impurities, no matter how weak, can easily pin down them in space at 
sufficiently low temperatures, truly breaking the translational invariance, and therefore causing the diminishing of the 
degeneracy (and thus the entropy). The thermopower should then quickly deviate its high-temperature saturation 
value and vanish as T ^ 0. In this regime, the Mott VRH will dominate the charge transport, as to be discussed 
below. 



B. Resistivity 



1. Mott variable range hopping at low temperatures 



Experimentally, the cuprate superconductors have universally exhibited the localization of charge carriers at low 
temperatures, in the low-doping regime of S < 0.05. The resistivity can be well fit null 13 by the following Mott 
VRH formula 

PM - e(-) (30) 

with 7 ~ 3—4 and Tq ^ 10^ K at T ^ 0, usually applicable to a doped semiconductor. This implies a strong localization 
of the doped holes in this regime. But lightly doped cuprates by no means resemble a doped semiconductor. The 
strong Coulomb interaction makes it a doped Mott insulator, in which the doped holes interact strongly with the spin 
background. As the result, they can be self-trapped at low doping even without any disorder as described before. 

So the self-localization of the doped holes in the lightly doped Mott insulator will provide an intrinsic mechanism 
to explain the localization phenomenon generally observed in the cuprates. Disorder or impurities, on the other 
hand, should only play a minor role in such a system. As noted above, in the presence of disorder, the spontaneous 
translational symmetry breaking of the lightly doped Mott insulator (with a lot of degeneracies) can easily become 
truly translational breaking, as the self-localized holes, without the penalty from the kinetic energy, are energetically 
in favor of staying near the impurities. Therefore, the low-doping phase can also be regarded as a strong Anderson 
localization system at low temperatures, even though the presence of disorder or impurities may not be really strong. 
In other words, the impurity effect will get "amplified" by the Mott physics at low doping. 

Recall that a holon has its own bare hopping term, governed by Hh in Eq. , which in the continuum limit reduces 

to 

with an effective mass = 2 f^nd A* as the continuum version of the gauge field . So the holon is expected to 
hop around based on Hh and is bound to the induced antimeron by the attractive potential H16|l . At low temperatures, 
A* may be neglected as spinous are in RVB pairing ,15J . Then the Schrodinger equation for a hole-dipole can be 
written down by 

+ Vip = E^i^. (32) 

2m/i 
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Define V'(r) = ipif^ 4')- To compute the radial component of the wave function we note that asymptotically at large r 
the Schrodinger equation reduces to a form whose (radial) solutions can be expressed as 



e (33) 

where oq — 



Then one may estimate the transition probability of the holon between any two adjacent antimerons ("impuri- 
ties"), located at i and j, based on Ea. H33|l . It is given by 



r.-exp(-^AL^-^) (34) 

where is the distance and is the on-site energy difference between two "impurity" sites. Except for the factor 

^In this formula is essentially the same as in the original Mott theory. As the temperature is lowered, the motion 

between neighboring sites becomes more difficult due to the lack of appropriate energy differences. Consequently, 
it is more likel y fo r the carriers to hop to a more distant site if this means that the energy difference is less. It is 
known |2a. l29l l3(]| as the Mott VRH. Except for a logarithmic correction, the resistivity for the 2D Mott VRH can 
be determined according to Eq. (|34|l by the following expression 



p{T) - e^ — ) V (35) 

13.8 



where Tq is a characteristic temperature given by Tq ~ tt^ ■ Here Dq is the energy density of the impurity states 



which is assumed to be constant in the VRH regime. 

It has been well known that once the anisotropic 3D is considered, the exponent 7 in the Mott VRH conductivity 
generally will be changed from 3 to 4, with a modified IqEI- Furthermore, the interacting effect between the holons, 
given by V12 = —ttJs In ili^^lsl^ j^ag been ignored here, which can also modify the exponent in the VRH theory. 



2. Crossover to deconfinement at high-temperatures 

The "deconfining" temperature Tdc will represent a characteristic temperature beyond which holons are deconfined 
from the antimerons. Once the holons are unbound from their antimeron partners and move freely at T > Tdc, 
their transport will be solely governed by the hopping term of the phase-string model H31|l . Here the interaction 
between holons and antimerons becomes irrelevant as F — > in the above RG analysis. Instead we must consider 
the contribution from the gauge fluctuations of A'* in Ea. (|31|l . which will play an essential role for scattering at high 
temperatures. 

Note that A* satisfies the following condition: 



dr . A'^ = ±7r ^ (nf^ - nf^) ^ $c. (36) 
We can estimate the strength of the fluctuations of A'* by defining 



ti7 = ^(A$^), (37) 

where $n denotes the flux per plaquette: 



$n = ±2^i^5f (38) 



4 

□ 



according to Ea. l|36|l . At very high-temperature limit, one may neglect the NN spin-spin correlations such that 

2^^< (5f)2 > 



= T:VT~Sr^TT, (39) 



which implies very strong flux fluctuations per plaquette in the high-T limit. 
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FIG. 4: The dielectric constant k as a function of doping calculated at different temperatures, k = 1 at half-filling, and the 
dashed lines indicate the linear dependence of the doping concentration. In the main panel, n diverges at Sdc ~ 0.038 and in 
the inset Jde = 0.008. 

As a matter of fact, if the fluctuations of A'* is treated in the quasistatic limit with an annealed average over static 
flux distributions in Ea. l|31|) . the transport properties are the same as those studied in Ref. 32]. In particular, the 
scattering rate has been found 3^] 

- ^ 2fcBT, 

T 

if > 7r/2, which is satisfied in our case according to Ea. (|39|l in the high-T limit. Corresponding, the resistivity is 

p - T at T > 1. (40) 

Namely, the charge transport in high-T "deconfining phase" will generally follow a strange-metal behavior due to the 
scattering between the holons and gauge field A'' in Ea. l|31|l . 

In the crossover from the low-T VRH behavior (|35|l to the high-T linear temperature behavior (|40|l . one expects to 
see a minimal resistivity pmin- We point out that generally the "deconfining temperature" Tdc does not necessarily 
coincide with the characteristic temperature of Pmin- The latter may occur at a lower temperature as the fluctuations 
in A" can already become important when the Neel temperature T/v is approached from below. In the above discussion 
of the VRH resistivity, such a scattering effect has been neglected at low-T, which should lead to the enhancement 
of the resistivity once it becomes important. Another possibility is that the holon induced antimerons may start to 
move with the increase of temperature, resembling the "flux-flow" in a superconducting phase, which also may lead 
to a qualitative change of the resistivity. These possibilities in the intermediate temperature regime are beyond the 
scope of the present work. 

C. Dielectric constant 

Another interesting prediction of the self- localization of doped holes is the existence of a large and doping- dependent 
dielectric constant as each hole is a dipolar object of a bound state of a holon and a localized antimeron. One thus 
expects that the dielectric constant increases linearly with doping initially and flnally diverges as the deconfining 
point is approached with increasing doping concentration. 

The dielectric constant defined in the RG analysis in Sec. II C can be written by k{1) — X{l)(3q^ /2-k. We can 
determine k = k{1 = oo) = X{1 = oo)(3q^ /2tt by Eq. H24|) numerically (taking k{1 = 0) = 1). Fig. 4 and the inset 
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show the results for k— latr = 0.1|^ and T — respectively. At half-filling, k = 1 as no contribution from the 
hole-dipoles. We see that with the increase of the hole concentration, the dielectric constant grows linearly with S 
at first, then deviates the linearity shown by the dashed line, indicating the increase of the dipole size. Eventually it 
diverges at a critical doping concentration 6dc — SdciT), beyond which the dipoles will collapse and free holons will 
be released. 

The general trend of the calculated in-plane dielectric constant, shown in Fig. 4, qualitatively agrees with 
the experimental measurements I5I, !^ in the low-doping cuprates. Indeed, for the lightly-doped cuprates, a large, 
doping-dependent dielectric constant has been observed in the a6-planc, which increases with the hole concentration 
S, initially linearly then becoming divergent at some higher concentrations 5, 6J. In contrast, the out-plane (c-axis) 
dielectric constant shows no essential change as a function of 6. 

In a doped-semiconductor picture, a dielectric constant contributed by the doped holes can only be significant when 
the holes are bound to impurities, which is in the localized regime at low temperatures. But such a dielectric constant 
should be usually anisotropic 3D-like rather than pure 2D-like as the experiment revealed. Especially it is difficult to 
explain why the dielectric constant should diverge in the a6-plane while remains constant in the c-axis. Furthermore, 
if the majority of the holes remains bound to impurities, it is hard to reconcile with the large saturated thermopower 
observed at T > 100 K. 

In contrast, the dipolar structure of the doped holes in the present framework can naturally lead to a large dielectric 
constant in the ab-plane, no matter whether the hole-dipoles are pinned down by impurities or not, as long as the 
hole dipole composites remain stable. The hole dipoles are presumably de-pinned from impurities at T > 100 K in 
our theory, since a large thermopower has been seen in experiment. With further increasing temperature, the dipolar 
structure will eventually collapse and the holon will become deconfined from its antimeron partner. Consequently the 
large dielectric constant and thermopower should be both quickly reduced above Tde, where the resistivity also starts 
to behave like a strange-metal as T becomes sufficiently high. 



IV. CONCLUSIONS 



In this paper, we studied the motion of doped holes in a spin ordered background at low doping. Based on the phase- 
string model, we demonstrated that the holes will get self-localized in space, leading to a spontaneous translational 
symmetry breaking without the presence of disorder or impurities. 

This novel property is an important consequence of the Mott insulator at low doping, described by the phase- 
string model as the low-energy effective description of the t — J Hamiltonian. The doping effect and the interplay 
between charge and spin degrees of freedom are characterized by a unique gauge structure with a mutual duality. 
At low doping, the spinon condensation forces a "confinement" on the holons, making the latter self-localized and 
resulting an insulator with an AFLRO or spin glass. This is in contrast to the higher-doping phase, where the holon 
condensation forces a "confinement" on the spinon part, resulting in a superconducting phase coherence [T^ . 

We found strong experimental implications based on the self-localization of holes. Large and doping-dependent 
thermopower can naturally explain the experimental data which had been very hard to understand by conventional 
theories. A large and doping-dependent in-plane dielectric constant indicates a composite structure of the holes and 
provides a unique explanation of the experimental observations, which otherwise are very difficult to comprehend. 
Furthermore, the low-T VRH resistivity observed experimentally was interpreted as the direct consequence of the 
self-localization with disorders playing a minor role, which explains why the critical doping of the delocalization in 
the cuprates is universally around 5c ~ 0.05 at T = 0, not very sensitive to the density of disorders in the samples. 
The phase-string model also naturally shows how the resistivity evolves into a strange-metallic linear-T behavior at 
sufficiently high temperatures above the delocalization temperature. Most importantly, we wish to emphasize that all 
these peculiar experimental properties were shown to be consistently explained within a single theoretical framework. 

It should be noted that many results in this paper arc only correct for a homogeneous phase. There may exist 
another possibility, namely, the stripe instability 10] in the phase string model, which can result in an inhomogeneous 
phase. Since holes are self-localized, their kinetic energies are suppressed such that the potential energy will become 
predominant. The dipole-dipole interaction might cause stripe instability at low temperatures, with hole-dipoles 
collapsing into a one-dimensional line-up (stripe) 'itf]. The pinning effect of disorders may stabile the homogeneous 
phase at low-T, so does the long-range Coulomb interaction. But it would be very interesting to incorporate the 
inhomogeneous tendency in various dynamic properties at low doping in future investigations. 
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